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We have addressed the microscopic transport mechanism at the switching or “on-ofF” transition in 
transition metal dichalcogenide (TMDC) field-effect transistors (FET), which has been a controver¬ 
sial topic in TMDC electronics, especially at room temperature. With simnltaneous measurement 
of channel condnctivity and its slow time-dependent fluctuation (or noise) in ultra-thin WSe2 and 
M0S2 FETs on insulating Si 02 substrates, where noise arises from McWhorter-type carrier number 
fluctuations, we establish that the switching in conventional backgated TMDC FETs is a classical 
percolation transition in a medinm of inhomogeneous carrier density distribution. From the ex¬ 
perimentally observed exponents in the scaling of noise magnitude with conductivity, we observe 
unambiguous signatures of percolation in random resistor network, particularly in WSe2 FETs close 
to switching, which crosses over to continuum percolation at a higher doping level. We demonstrate 
a powerful experimental probe to the microscopic nature of near-threshold electrical transport in 
TMDC FETs, irrespective of the material detail, device geometry or carrier mobility, which can be 
extended to other classes of 2 D material-based devices as well. 


The expanding family of atomically thin layers of semi¬ 
conducting TMDG materials for electronic [1-7], opto¬ 
electronic [8-14], valleytronic [15-19] and even piezo- 
electronic [20] applications, has defied a generic frame¬ 
work of electron transport because of diverse material 
quality, channel thickness dependent band structure, 
dielectric environment, and nature of substrates. A 
wide variety of physical phenomena ranging from vari¬ 
able range hopping [5], metal-insulator transition [7] to 
classical percolative charge flow through inhomogeneous 
medium [21,22] were reported for M 0 S 2 , the implications 
of which often lead to a conflicting microscopic scenario. 
At low carrier density, for example, hopping via single 
particle states trapped at short-range background poten¬ 
tial fluctuations (~few lattice constants [23]) is incompat¬ 
ible to the observation of classical percolative conduction 
that requires long-range inhomogeneity in charge distri¬ 
bution, created when linear screening of the underlying 
charge disorder breaks down [24-27]. Observations of 
metal-insulator transition [7,22] have added to this con¬ 
troversy with both many-body Coulomb interaction as 
well as classical percolation transition being argued as 
possible driving mechanisms. With the emergence of 
new TMDC-based FETs, in particular WSe 2 [19,28-35], 
WS 2 [36-39], and MoTe 2 [6], the importance of a generic 
probe to the microscopic details of density (or gate volt¬ 
age) dependent charge distribution that affect the mobil¬ 
ity, subthreshold slope and other performance markers in 
TMDC electronics, is paramount. 

Unlike variable range hopping in the strongly local¬ 
ized regime, however, identifying charge percolation in 
FET channels is less straightforward. Theoretically, clas- 
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sical percolation-limited conductivity in metal-insulator 
composites is characterized by established critical expo¬ 
nents [21,24,25,40-44], but the experimental difficulty lies 
in accurately determining the fraction p of the conduct¬ 
ing region, or “puddles”, embedded inside the insulat¬ 
ing matrix. Hence despite compelling evidence of long 
range inhomogeneity in the charge distribution in M 0 S 2 
FETs [21,22], its manifestation in transport remains in¬ 
direct and confined only to low temperatures. A way 
to circumvent this difficulty involves measuring the low- 
frequency noise, or l//-noise, in the channel conductivity 
CT, which also scales with p with independent characteris¬ 
tic scaling exponents [40,42], and diverges at the perco¬ 
lation threshold {pc) [41]. A direct relation between the 
normalized noise magnitude and cr thus eliminates 
the necessity to know either p or pc accurately. In such 
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where the scaling exponent v assumes universal values 
depending on the nature of percolation. In the 2D lattice 
models of random resistor network, for example, v « 0.86 
or 1.5, [REF: [40,45]] whereas in the continuum perco¬ 
lation framework, such as the “swisscheese model” [44], 
v fs 3.2 for random void (RV, insulating voids in a con¬ 
ducting background) and « 0.87 for inverted random 
void (IRV, weakly connected conducting regions in an 
insulating matrix), respectively [41,42]. In this work we 
have investigated the switching process in TMDG FETs 
by simultaneously measuring Ng- and a as the electri¬ 
cal transport is tuned from the strongly insulating to the 
quasi-metallic regime using a gate voltage. While noise in 
all devices was found to primarily originate from carrier 
number fluctuations at the channel-substrate interface, 
the key result is the observation of an unambiguous scal¬ 
ing of Nc with (j, suggesting classical percolation in spa¬ 
tially inhomogeneous medium at the onset of conduction. 
Similar behavior for both WSe 2 and M 0 S 2 FETs imply 
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TABLE I: Channel Length (L), Channel Width (W), details of electrical contact and layer number for the WSe 2 FET devices 


Device Name 

L (pm) 

W (pm) 

Contact Details 

No. of layers 

WSeML_2 

2.8 

3 

Ti/Au 10nm/50nm 

5 

WSeBL 

1.5 

5.3 

Ti/Au 10nm/50nm 

2 

WSeML_3 

1.8 

3 

Ti/Au 10nm/50nm 

5 

MoSSL 

3 

4 

Cr/Au 10nm/50nm 

1 

MoSSL.l 

0.876 

0.861 

Cr/Au 10nm/50nm 

1 

MoSSL_2 

0.562 

1.023 

Cr/Au 10nm/50nm 

1 


that percolative switching could be generic to backgated 
TMDC FETs irrespective of material and device-related 
parameters. 

We chose WSe 2 FETs as the primary experimental 
platform due to the following reasons: First, WSe 2 is 
an emerging TMDC FET material with several desirable 
properties ranging from high carrier mobility due to low 
effective mass of the carriers, ambipolar conduction and 
superior chemical stability compared to sulphides [12,13, 
28,32,33,35,37]. Second, in spite of the progress in stan¬ 
dard electrical transport properties [19,28-35], the origin 
and magnitude of intrinsic 1//-noise, a crucial perfor¬ 
mance limiting factor in electronic device applications, 
in WSe 2 FETs is not known, and third, given the recent 
studies of noise in M 0 S 2 FETs [46-50], similar studies in 
WSe 2 allows identification of the generic aspects of noise 
processes in TMDC FETs, which in turn, provides cru¬ 
cial insight to microscopic details of charge distribution 
and disorder. 

The experiments were carried out in FETs fabricated 
from ultra-thin films of WSe 2 (and M 0 S 2 ) obtained 
via mechanical exfoliation on p^'^-Si/(285 nm)Si02 
substrates in the conventional backgated geometry 
(schematic and typical device image in Fig. la) (De¬ 
tails of the measurement technique provided in Supple¬ 
mentary Information). The details of the devices can 
be found in Table I. In order to verify that the struc¬ 
tural integrity of the WSe 2 channels were unaffected by 
the device fabrication process, we carried out Raman 
spectroscopy on, (1) original block of WSe 2 which was 
subsequently mechanically exfoliated to form the device 
WSeML_2, (2) channel region prior to device fabrica¬ 
tion, and (3) the same after device fabrication. The nor¬ 
malized Raman spectra (Fig. lb) remains essentially un¬ 
changed, with characteristic Raman modes for multilayer 
WSe 2 including those at « 247.3 cm“^ 249.5 cm~^ 

(A'l), and « 308.5 crn’i {AlJA'l) [51], implying no 
detectable structural damage due to the exfoliation or 
the device fabrication processes. Similar characteriza¬ 
tion for the M 0 S 2 channels and detailed surface char¬ 
acterization with atomic force microscopy can be found 
in the supplementary material. Measurement of both 
standard transport and noise were performed in a two- 
probe geometry via ac lock-in amplifier technique, where 
the latter was recorded by digitizing the amplifier out¬ 


put by a analog-to-digital converter over a bandwidth of 
« 0.01-10 Hz [47]. 

In TMDC FETs, thermal annealing and environmen¬ 
tal exposure have drastic effects on electrical mobility, 
threshold voltage and other parameters [35,37,52]. Re¬ 
cently, we have shown [47] that thermal annealing can 
also reduce the noise magnitude in M 0 S 2 FETs by over 
two orders of magnitude, through improved transparency 
of the electrical contacts [37,52]. We found a strong ef¬ 
fect of annealing on the electrical transport in our WSe 2 - 
FETs. In Fig. Ic-f, we illustrate the effect of annealing on 
both transport and noise characteristics in WSe 2 FETs, 
as well as the sensitivity of these parameters to envi¬ 
ronmental contamination. The inset of Fig. Ic shows the 
dependence of the source-drain current {Isd) on gate volt¬ 
age Vg in device WSeML_2 under three conditions: (a) in 
air/atmosphere prior to annealing (labelled as Air), (b) 
in vacuum (rj 10~® mbar) prior to annealing (labelled as 
Vacuum), and (c) after one hour of annealing at 400 K 
(in vacuum) (labelled as Vac. Anneal). The annealing 
conditions were kept the same for all WSe 2 and M 0 S 2 de¬ 
vices reported here. Apart from the evident increase in 
Isd on annealing, typically to about ~ 100 nA at largest 
experimental Vg (Fig. Ic) {Vsd = 100 mV), we also note: 
(i) First, a shift of the turn-on voltage Vqn, indicated by 
the dashed vertical lines, towards negative values. This 
has been observed earlier in WS 2 [REF: [37]] and M 0 S 2 
[REF: [53]] channels, and implies desorption of gas and 
water dopants. (Here, Vqn is calculated by extrapolat¬ 
ing the \/T^ vs Vg plot in the ON state of the transis¬ 
tor and finding its intercept on the Vg axis [54].) This 
makes p-type conduction in our WSe 2 channel inaccessi¬ 
ble within the experimental range of Vg, although it al¬ 
lows a wide operating range in the n-doped regime, (ii) 
Second, nearly four decades of increase in the field-effect 
mobility (p = {llCox)d(jldVg, where Cox is the gate ca¬ 
pacitance per unit area) (Fig. Id) to /x « 0.5-5 cm^/V.s, 
in agreement with previous reports [37]. (iii) Third, lin¬ 
ear and symmetric Igd ~ Vsd characteristics at low source- 
drain voltages (up to \Vsd\ S 100 mV) (supplementary 
material), confirming the ohmic nature of the contacts. 
All electrical measurements were limited to the linear 
transport regime, where the fluctuations in Igd reflects 
that of cr. Thermal annealing in vacuum had a dra¬ 
matic effect on the intrinsic low-frequency l//-noise mag¬ 
nitude in WSe 2 FETs, particularly at large values of Vg 
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FIG. 1: Device structure and transport characterization, (a) Device Schematic and optical micrograph of a typical WSeg 
FET. (b) Raman Spectra for bulk WSeg, WSeg flakes after exfoliation (Pre-fabrication) and following device processing (Post¬ 
fabrication). (c) Transfer charateristics of three WSeg FETs after annealing. The dashed lines indicate the turn-on voltage 
(see text). (Inset) A comparison of the transfer characterisitcs in air (labelled ‘Air’), in vacuum before annealing (labelled 
‘Vacuum’) and in vacuum after annealing (labelled ‘Vac. Anneal’) for the device WSeML_2. (d) Comparison of the field effect 
mobility before (BA) and after (AA) annealing curves for three devices WSeML_2, WSeML_3 and WSeBL. (e) Normalized 
power spectral density (PSD) of current fluctuation showing a l//“-type frequency dependence. Lower inset: Vg-dependence of 
frequency exponent a both before and after annealing. A large reduction in the normalized time dependent current fluctuations 
is observed after annealing (upper inset), (f) Comparison of normalized noise magnitude for WSeML_2 before and after 

annealing as a function of gate voltage Vg measured from the turn-on voltage Vbiv- Inset: Area-scaled normalized noise 
magnitude after annealing for two devices WSeML_2 and WSeML_3. 


(high doping). This is illustrated in the upper inset of 
Fig. le, which shows the normalized time-dependent cur¬ 
rent fluctuations Slsdlhd at the same gate bias before 
and after the annealing process. The overall shape of the 
power spectral density £'/(/) of the fluctuations, which 


depends on frequency / as Si{f) oc l//“ with a rs 1 
(Fig. le), is insensitive to the annealing state or Vg (lower 
inset of Fig. le). The main effect of annealing, how¬ 
ever, lies in the V),-dependence of the normalized noise 
magnitude and shown in Fig. If with the data 
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FIG. 2: Integrated noise and number fluctuation. Plots of integrated power spectral density (markers), calculated McWhorter 
noise Nmcw (olive green colored line) using substrate trap density Da as fitting parameter (Equation 2) and source drain 
current (navy blue line) for (a) WSeBL, (b) WSeML_2, (c) WSeML_3. The remarkable agreement between the calculated 
McWhorter noise with the experimental noise magnitude indicates channel carrier number fluctuations as the source of flicker 
noise in our WSe 2 channels. 


from device WSeML_2. The pre-anneal V^-dependence 
of is weak with a monotonic decrease as Vg is 

increased, whereas the variation in becomes non¬ 

monotonic after annealing, and can be smaller than the 
pre-anneal noise by more than two orders of magnitude at 
large Vg [47,55]. Quantitatively, the total measured noise 
can be decomposed into noise from the contact (nor¬ 
malized noise Scon) and the channel (normalized noise 
Sch) regions, so that Sj/l'^^ = ij'^Scon + (1 - with 

V = GchKGch + Gcon), where Gch and Gcon are the chan¬ 
nel and contact conductances, respectively [46]. Prior to 
annealing, rj ^ 1, and the contact noise is expected to 
dominate, while the post-annealing (77 < 10 “'^) noise is 
determined by the channel region. To confirm this we 
examined the dependence of the post-anneal normalized 
noise magnitude on the channel area (A) of two devices of 
different sizes from the same WSe 2 flake. As shown in the 
inset of Fig. If, from the devices collapse at all 

Vg, implying oc If A, which is expected when noise 

originates from the channel region [56]. (The 1/A depen¬ 
dence was not observed in the pre-annealed noise.) In 
the remainder of this paper we have confined our discus¬ 
sion only to the post-anneal noise behavior in the TMDC 
channels. Moreover, in order to analyze the noise data 
without referring to any specific frequency, we have com¬ 
puted the “variance” or the integrated power spectral 
density as = {6a'^)/{af = {l/{Isdf) f Si{f)df, where 
the integration was carried out within the experimental 
bandwidth. 

In order to explore the microscopic mechanism of 
noise in WSe 2 FETs, we note the strong peak close 
to the steepest rise in the Igd - Vg transfer character¬ 
istics at a device-dependent gate voltage Vg in all WSe 2 
FETs (Fig. 2). This naturally suggests the noise from 
McWhorter carrier number fluctuations [57], where the 
fluctuation in I^d is due to trapping and detrapping of 


free carriers at the channel-oxide interface. This is equiv¬ 
alently represented by the flat-band gate voltage fluc¬ 
tuations, and thus resulting in current noise being pro¬ 
portional to the transconductance gm = dIsd/dVg of the 
channel. To establish this quantitatively, the noise mag¬ 
nitude can be written as [46-50,57], 


Nmcw 


{Isd? 


eAe'^ksT (gm]^ 

agi,k ■ 


( 2 ) 


where Da ~ 10 ^® - 10 ^° eV"^cm“® is the density of trap 
states in Si 02 , and k « 10 ® m“^ is the inverse of elec¬ 
tronic wave function decay scale inside the oxide. Using 
Da as the fitting parameter, and experimentally deter¬ 
mined gm, we have compared the Ug-dependence of the 
observed noise magnitude and that expected from 
the McWhorter model (Equation 2) for all WSe 2 EETs 
(Fig. 2). The excellent agreement with realistic range 
of values of Da for Si 02 establishes the carrier number 
fluctuations at the channel-Si 02 interface as the predom¬ 
inant source of noise in our WSe 2 FETs. 

Subsequently, we have examined the dependence of 
on the channel conductivity a, which shows very sim¬ 
ilar trends for all three WSe 2 devices. In Eig. 3a, c 
and d, three distinct regions can be clearly identified: 
(i) Region I: The low-conductivity regime at Vg < Vg, 
where increases from a finite value with increasing a 
(dashed lines) to the maximum at Vg which corresponds 
to CT 3 - 5 X 10~® S/n, (ii) Region II: &i Vg > Vg ex¬ 
tending over a ~ (0.5 3 ) x lO'"^ S/n (indicated as the 

shaded area in Fig. 3a-d) where oc cr^G-^^o. 2 ) ^ 

(iii) Region III at Vg » Vg and cr > 3 x lO”"^ S/n, where 
decreases rapidly as oc The three 

regimes are characterized by distinct T-dependences of 
cr, as illustrated for WSeML_2 in Fig. 3b, although other 
devices show very similar behavior (supplementary ma- 
















5 


a 


c 



0.004 0.006 0.008 0.010 





N„ PC (RV) 


0 

’ 1 

1 • 



Fermi sea 


] w 

N„oc a° 

‘*HiRV)/o-^HRR) 

^_J 

■ » 

Void 





FIG. 3: Noise - conductivity scaling and percolation. Dependence of noise magnitude on conductivity cr and the gate 
voltage dependence of a (insets) for (a)WSeML_2, (c)WSeML_3, (d)WSeBL. Different noise scaling regimes are indicated by 
solid lines with respective exponents {v). The shaded region outlines the transition region with u = 1.5. (b) Temperature 
dependence of a of WSeML_2 indicating transition from localized to almost temperature independent transport, (e) Schematic 
of the “island-and-sea” representation of charge distribution for different fermi level positions as a function of gate voltage. 


terial) [35,37,52]. The weak T-dependence of a in Re¬ 
gion III signifies quasi-metallic electron transport. In 
Region I and Region II, i.e. around the switching tran¬ 
sition (cr < 10~^ - I0~® S/n), the insulating behavior 
becomes progressively stronger with decreasing V),, al¬ 
though the signature of variable range hopping devel¬ 
ops only at much lower Vg i.e. stronger localization 
(see Fig. 3b and supplementary material for more de¬ 
tail), which makes identifying the transport mechanism 
at switching from T-dependence of a alone a very difficult 
task. 

The exponents v in N„ - a scaling, however, provide 
a direct evidence of a percolative transport, with the 
noise peak signifying the percolation threshold a.t Vg 
Vg [40,42,58]. Below the percolation threshold {Vg < Vg-. 
Region I), called the “dielectric regime”, transport occurs 
by hopping or tunneling through disconnected metallic 
puddles [59], and approaches a finite device-dependent 
magnitude at low Vg away from the threshold. In the 
critical (metallic) regime immediate to the percolation 
threshold {Vg > Vg-. Region II), v ps 1.5 corresponds to 
percolation in random resistor network (RR in Fig. 3(e)) 
with a broad distribution of intersite resistance [40]. On 
further increase in Vg, the crossover to iz « 3 in Region III 
marks the onset of continuum percolation with random 
voids, signifying transport in delocalized Fermi sea in the 
presence of insulating islands [42]. 

Similar variation of v, and thereby the nature of per¬ 
colative transport, has been observed earlier in a num¬ 
ber of metal-insulator composites with varying fraction 


of the metallic component (see, for example, REF [58]). 
In our case, it allows a unique route to monitor the evo¬ 
lution of carrier distribution with gate voltage, or Fermi 
energy, especially at switching. Switching by percola¬ 
tion transition within the “island-and-sea” description 
in Fig. 3e, has been reported in several different classes 
of low-dimensional systems, including dilute 2D elec¬ 
tron or hole gases in high mobility semiconductor het¬ 
erostructures [24,25,60,61], silicon inversion layers [26], 
perovskite oxides [62], graphene nanoribbons [27,63], or 
even M 0 S 2 [21,22] and WS 2 [37] FET devices. How¬ 
ever the signatures of percolative transport in these sys¬ 
tems are usually indirect, and depend on the observa¬ 
tion of metal-insulator transition [24,25], or characteris¬ 
tic Isd - Vsd scaling [64,65], which can also be limited to 
low temperatures and/or suffer from nonequilibrium ef¬ 
fects. The noise-conductivity scaling constitutes a new 
probe to percolative transport in 2D TMDC materials, 
where the divergence of noise associates the switching 
transition of the FET to the percolation threshold. In 
order to verify the generality of the observed Na - a 
scaling in other classes of TMDC FETs, we have car¬ 
ried out noise experiment in similarly prepared back- 
gated single layer M 0 S 2 FETs on Si02 substrate (see 
Table I for details). Three different devices with field- 
effect mobility upto p ~ 20 cm^/V.s (inset of Fig. 4a) 
were measured across the switching transition covering 
nearly five decades of change in conductivity (Fig. 4a). 
Fig. 4b (inset) illustrates the typical gate voltage depen¬ 
dence of iVcr in one of the devices (MoSSL_l), where we 
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FIG. 4: Transport and Percolation Noise in M 0 S 2 FET. (a) Gate voltage dependence of conductivity (a) and field effect mobility 
(fi) (inset) for three M 0 S 2 devices MoSSL, MoSSL_l and MoSSL_2. (b) Plot of normalised noise magnitude (N^-) vs gate bias 
(Vg) for MoSSL_l (inset). The dashed line indicates the noise contribution from the contacts. The main panel shows the gate 
voltage dependence of corrected for contact noise {Ncontact) and the calculated McWhorter noise (solid line) using trap 
density Du as a fitting parameter for device MoSSL_l. (c) Scaling of No - Ncontact with a for three M 0 S 2 devices showing 
characteristic percolative scaling exponents. The shaded region (in (a) and (c)) indicates the transition from the Inverted 
Random Void (u rj 0.85) to Random Void (u « 3) mode of transport. 


observe the noise magnitude to saturate at large doping 
{Vg > -40 V). We have shown previously [47] that such 
a saturation arises from remnant effects of the contacts, 
and contributes additively to the total measured noise. 
Subtracting the contact noise Ncontact (dashed line in 
Fig. 4b inset), obtained by fitting the V],-dependence of 
No citVg> -35 V, we recover the overall agreement be¬ 
tween the channel noise {N^ - Ncontact) and the calcu¬ 
lated McWhorter noise Nmcw (Equation 2) over nearly 
five decades suggesting that the noise in the annealed 
M 0 S 2 FETs has a very similar microscopic origin as the 
WSe 2 devices. 

Fig. 4c shows the scaling of - Ncontact with the 
channel conductivity cr for all three devices (shifted ver¬ 
tically with respect to each other for clarity). The scaling 
behavior is very similar in all cases where i/ rs 0.85 ± 0.15 
at low a and crosses over to the continuum 2D perco¬ 
lation with random void (z/ « 3) above a characteristic 
conductivity Uc ~ 5 x 10~^ S/n. The McWhorter noise, 
estimated from experimentally measured gm (shown in 
Fig. 4c for MoSSL_2 with open circles) behaves in identi¬ 
cal fashion over the entire range of a ensuring the consis¬ 
tency of the analysis. Interestingly, the magnitude of CTc 
is similar for WSe 2 and M 0 S 2 FETs in spite of different 
material and channel mobility. Importantly, M 0 S 2 chan¬ 
nels only show a divergence of noise with decreasing Vg, 
and no peak is observed even down to the smallest mea¬ 
sured (T (~ 10"^^ S/n). This suggests that unlike WSe 2 
FETs, conductivity below percolation threshold in M 0 S 2 


FETs is too small to be measurable, possibly due to a 
predominantly long range nature of disorder. Hence the 
observed ly (rs 0.85) at a < Uc in M 0 S 2 is more likely to 
be due to continuum percolation in inverted random void 
(IRV, Fig. 3e) configuration, rather than random resistor 
network. Nevertheless, the N^ - a scaling in M 0 S 2 FETs 
with universal exponents suggests that percolative trans¬ 
port at switching could be generic in TMDC FETs, al¬ 
though the nature of percolation itself (discrete random 
resistor or continuum), is likely to depend on material 
details, channel mobility, disorder landscape, and other 
factors. 

In conclusion, we have employed low-frequency I//- 
noise in electrical conductivity to explore the nature of 
charge transport in FET devices with ultra thin TMDC 
channels, close to the switching transition at room tem¬ 
perature. The noise in both WSe 2 and MoS 2 -based FETs 
arise from carrier number fluctuations, in agreement with 
previous experimental reports. Importantly, a unique 
scaling of the noise magnitude with channel conductiv¬ 
ity in WSe 2 FETs provides direct evidence of percolative 
electron flow at switching which, in WSe 2 EETs, crosses 
over from the random resistor network to continuum per¬ 
colation with increasing carrier density. Similar behavior 
in M 0 S 2 FETs suggest that percolative transport in in¬ 
homogeneous charge distribution at switching could be 
generic to other TMDC channels as well. 

We acknowledge the Department of Science and Tech¬ 
nology (DST) for a funded project. 
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TABLE SI: Ghannel Length (L), Channel Width {W), details of electrical contact and layer number for the 

unannealed WSe 2 FET devices 


Device Name 

L (/im) 

W (fim) 

Contact Details 

No. of layers 

WSeMLT 

0.5 

3.4 

Au 40 nm 

4 

WSeTL 

0.7 

0.5 

Ti/Au 10nm/50nm 

3 


I. EXPERIMENTAL DETAILS 

All devices were measured in a two probe geometry at a pressure of IE—6 mbar inside a cryostat. For transfer 
characteristics source drain bias to the device was supplied from LIASR830. The output current at the frequency of 
excitation was measured using the lock-in technique. Gate voltage was supplied using a Keithley 2400 DC sourcemeter. 
During I-V characteristics both gate voltage and source drain bias was applied using Keithley 2400 DC sourcemeters. 
For noise measurements, a bias voltage of 100 mV at a frequency of 220 Hz was applied to the sample. Time series 
data was obtained using NI DAG (NI USB-6210) at a rate of 1000 samples per second. The power spectrum was 
calculated by performing a fast fourier transform on this data. The Power spectral density follows 1// law of flicker 
noise Fig: 2{h). 



2 


II. BIPOLAR TRANSPORT IN WSea FETs 

One of the major advantages of tungsten based TMDCs WS 2 /WSe 2 over M 0 S 2 is the possibility of bipolar conduc¬ 
tion owing to their smaller bandgap and reduced fermi level pinning. Indeed, bipolar transport in WSe 2 FETs has 
led to the recent demonstration of electrostatically tunable integrated photovoltaic and photodetection devices [1-6], 
and also promises new designs of complementary digital logic by appropriate electrical contacting schemes [7]. In our 
experiments, some unannealed WSe 2 FETs (details in Table SI) did show bipolar conduction (Fig. SI) although the 
ON-state current in these devices were too small for any meaningful noise experiments. 



FIG. SI: (a) Transfer characteristics performed at 94.7 K for a range of source drain bias all showing bipolar channel 
conduction. The pinch off voltage for hole conduction was fairly high and gate oxide leakage prevented us from 
applying higher gate bias. Thus we see a small but prominent hole conduction. The mobility for electron and 
hole conduction was 0.6 and 0.1 cm^/(V.s) respectively, (b) Temperature dependent transfer characteristics showing 
bipolar transport at all temperatures from 110 K to 270 K. Data was obtained at a fixed source drain bias of 100 mV. 
This strong bipolar nature of conduction is a proof of the comparatively low Fermi level pinnning in WSe 2 when 
compared to M 0 S 2 . Maximum electron and hole mobilities for this device was 0.17 and 0.5 cm^/(V.s) respectively. 
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III. TEMPERATURE DEPENDENCE OF CONDUCTIVITY (a) AND VARIABLE RANGE HOPPING 

IN WSea FETs 


In Fig. S2, we have attempted to identify the microscopic transport mechanism, specifically the 2D variable range 
hopping or thermal activation. At low values of gate voltage (Vg) we observe an insulating behaviour which transforms 
to an almost temperature independent or slightly metallic behaviour at higher values of gate voltage. Neither two 
dimensional variable range hopping (Fig. S2) nor activated transport (shown for WSEML_2 in Fig. 3b) can describe 
the T-dependence of conductivity (a) adequately. However, at very low gate voltages {Vg < 30 V), Incr tends to 
become linear with (Fig- S2), probably indicating contributions from 2D variable range hopping. 


a 


b 


c 
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FIG. S2: (a),(b) and (c) Temperature dependence of conductivity for three separate WSe 2 FETs showing insulating 
behaviour for lower gate voltages which evolves into an almost temperature independent transport for higher gate 
voltages. The transition occurs through a percolative mechanism of transport as confirmed by the scaling of noise 
exponents as shown in the main body of the paper. 
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IV. EFFECT OF ANNEALING ON TRANSFER CHARACTERISTICS OF WSea FETs 


As stated in the main text, annealing leads to a significant increase in the drain current Igdi shifting of the turn-on 
voltage Von towards negative values and increase in mobility by almost four decades for WSe 2 FETs. These features 
were seen for all our annealed WSe 2 devices and is shown here for WSeBL (Fig S3(a)) and WSeML_3 (Fig S3(b)). 



FIG. S3: A comparison of the transfer characteristics in air (labelled Air), in vacuum before annealing (labelled 
Vacuum) and in vacuum after annealing (labelled Vac. Anneal) for the device WSeBL (a) and WSeML_3 (b). (c) 
Isd — Vsd characteristics of WSeML_2 showing linear dependence at low bias. 
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V. RAMAN SPECTROSCOPY AND AFM CHARACTERIZATION OF M 0 S 2 
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FIG. S4: (a) Raman characterization for M 0 S 2 flakes before and after device fabrication. Raman data for M 0 S 2 shows 
peaks for mode at « 385 cm“^ and Aig mode at « 405 cm“^. The raman peak at « 520 cm“^ is due to the 
silicon substrate. The position of the peaks match perfectly with previously reported data [8]. This leads us to believe 
that the crystal structure is not adversely affected by the fabrication process. We used a laser wavelength of 532 nm 
for our measurements, (b) Optical micrograph of an exfoliated M 0 S 2 flake. The white lines indicate the paths along 
which the height profile (panel (c)) for the M 0 S 2 and Si02 surfaces (shifted vertically for clarity) have been measured. 
The surface roughness in the channel region is « 0.604 nm whereas that on bare Si02 is « 0.65 nm. The similarity of 
these values indicate the conformity of the flakes to the Si02 surface. This percludes the presence of ruptures or other 
sharp structures in the channel region which might affect the transport properties adversely. Surface roughness was 
measured by computing the standard deviation of the height data in a given region of the channel or Si 02 surface. 
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VI. AFM CHARACTERIZATION OF WSez 



C 



FIG. S5: Optical micrograph (a) and a zoomed in AFM image (b) of a typical exfoliated WSe 2 flake. The height profile 
(c) (shifted vertically for clarity) has been measured along the paths (white lines) shown in panel (b). Roughness for 
both the flake and the Si02 surface was calculated to be « 1 nm. Hence, the WSe 2 flakes as in the case of M 0 S 2 ones 
(Fig S4) conform to the Si02 surface. 


[1] B. W. H. Baugher, H. O. H. Churchill, Y. Yang, and P. Jarillo-Herrero, Nat. Nanotechnol. 9, 262 (2014). 

[2] A. Pospischil, M. M. Furchi, and T. Mueller, Nat. Nanotechnol. 9, 257 (2014). 

[3] J. S. Ross, P. Klement, A. M. Jones, N. J. Ghimire, J. Yan, M. G., T. Taniguchi, K. Watanabe, K. Kitamura, W. Yao, 
D. H. Gobden, and X. Xu, Nat. Nanotechnol. 9, 268 (2014). 

[4] O. Lopez-Sanchez, D. Lembke, M. Kayci, A. Radenovic, and A. Kis, Nat. Nanotechnol. 8, 497 (2013). 

[5] D. J. Groenendijk, M. Buscema, G. A. Steele, S. Michaelis de Vasconcellos, R. Bratschitsch, H. S. J. van der Zant, and 
A. Gastellanos-Gomez, Nano Lett. 14, 5846 (2014). 

[6] M. M. Furchi, A. Pospischil, F. Libisch, J. Burgdrfer, and T. Mueller, Nano Lett. 14, 4785 (2014). 

[7] W. Liu, J. Kang, D. Sarkar, Y. Khatami, D. Jena, and K. Banerjee, Nano Lett. 13, 1983 (2013). 

[8] H. Li, Q. Zhang, G. C. R. Yap, B. K. Tay, T. H. T. Edwin, A. Olivier, and D. Baillargeat, Adv. Funct. Mater. 22, 1385 

( 2012 ). 









